Use of Vector Autoregressive model (VAR) to forecast hospital cases using multivariate data (with log transformation) with vaccination information in different age group and wastewater signal. 5 days of forecast.

Vector Autogressive model (VAR) is like ARIMA model which forecasts a time series based on the correlation between one timeseries that we are trying to forecast and its own lagged values. However, VAR being a multivariate model also allows use of multiple variables to be involved whose lagged values can also be used to forecast. However, unlike ARIMAX model which allows use of external regressors in the model in a one way causal relationship where the external variables can be used as regressors to forecast one univariate timeseries, VAR model allows two way causal relationship where the lagged values of one timeseries can cause changes in another timeseries and the lagged values of the latter timeseries can also cause changes in the former one. VAR model regresses each variable against its own lagged values as well as the lagged values of other variables in the dataset.

VAR model just like ARIMA also requires stationary time series to work with. In this notebook, we check for those assumptions or make the time series data stationary. The goal was to forecast hospitalization using the lagged values of other variables that might impact the hospitalization from vaccination data and the waste water signal. The most appropriate lag value was selected based on lowest AIC error metric from a loop of different lag values used.

The data has been log transformed and then divided into train and test set and the RMSE score, MAPE error metric has been evaluated.

Granger causality test has been performed to confirm whether there is a causal relationship between different time series in our data that is statistically signicant.

Initial Data Cleaning

Vaccine data, Cases data from the hospital, and the waste water signal data has been loaded, cleaned, and then merged into one final dataframe, final_data.

The variables have been log transformed.

Creating a time series dataframe from the original data using frequency of 365 for days.

Time Series Plots for those with 1 dose of vaccination in the age groups (5-40 years old) dataset below:

Time Series Plots for those with 1 dose of vaccination in the age groups (50 - 80 years old) and those with 2nd dose in the age groups (5 - 12 years old) in the dataset below:

Time Series Plots for those with 2nd dose of vaccination in the age groups (18 to 80 years old) and the plot for hospitalization and wastewater signal (N1_N2_avg).

To fit a VAR model on the data, the timeseries must be stationary with stable mean and variance over time. The time series above show definite trend and the mean is not constant over time. In wastewater and hospitalization, the variance do not seem to be constant over time. Therefore, they do not appear stationary. However, a statistical test to check whether the time series is stationary or not, Augmented Dickey Fuller Test has been performed.

Null Hypothesis H0: The time series is not stationary. Alternate Hypothesis H1: The time series is stationary.

From the p-values obtained below is evident that we can not reject the null hypothesis and our time series are not stationary.

#> $Percent_of_Ottawa_residents_5_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -1.3785, Lag order = 7, p-value = 0.8396
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_12_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -1.1089, Lag order = 7, p-value = 0.9204
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_18_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = 0.055666, Lag order = 7, p-value = 0.99
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_30_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -0.67072, Lag order = 7, p-value = 0.9728
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_40_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -1.0632, Lag order = 7, p-value = 0.9278
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_50_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -0.58684, Lag order = 7, p-value = 0.9777
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_60_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -0.79866, Lag order = 7, p-value = 0.9613
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_70_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -1.7984, Lag order = 7, p-value = 0.6623
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_80_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -2.8091, Lag order = 7, p-value = 0.2355
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_5_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = 1.0097, Lag order = 7, p-value = 0.99
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_12_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -1.7387, Lag order = 7, p-value = 0.6875
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_18_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -2.6859, Lag order = 7, p-value = 0.2875
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_30_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -0.648, Lag order = 7, p-value = 0.9748
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_40_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -0.59252, Lag order = 7, p-value = 0.9774
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_50_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -0.71766, Lag order = 7, p-value = 0.9686
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_60_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -0.76492, Lag order = 7, p-value = 0.9643
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_70_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -0.98347, Lag order = 7, p-value = 0.9407
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_80_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -1.8037, Lag order = 7, p-value = 0.6601
#> alternative hypothesis: stationary
#> 
#> 
#> $observed_census_ICU_p_acute_care
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -1.3428, Lag order = 7, p-value = 0.8547
#> alternative hypothesis: stationary
#> 
#> 
#> $N1_N2_avg
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -2.8537, Lag order = 7, p-value = 0.2167
#> alternative hypothesis: stationary

To make the time series stationary for fitting VAR Model, the first difference of the time series data has been taken by subtracting the current row from its previous row. Before differencing the time series the dataset has been divided into train and test datasets. The last days of the timeseries is part of the test dataset. The MTS package has been used for taking the differnce of the time series.

After taking the difference of the time series, the Augmented Dickey Fuller test is again performed on it to check if time series is now stationary. The p-value is statistically significant below 0.05. Therefore, we can say that the time series are now stationary.

#> $Percent_of_Ottawa_residents_5_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -6.004, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_12_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -3.812, Lag order = 7, p-value = 0.01862
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_18_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -5.2073, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_30_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -5.3235, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_40_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -4.7868, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_50_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -4.324, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_60_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -4.4641, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_70_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -4.4885, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_80_1_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -6.1205, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_5_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -3.5652, Lag order = 7, p-value = 0.03647
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_12_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -4.122, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_18_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -2.0296, Lag order = 7, p-value = 0.5647
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_30_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -5.3596, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_40_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -5.11, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_50_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -5.2996, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_60_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -5.0694, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_70_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -4.9772, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $Percent_of_Ottawa_residents_80_2_dose
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -5.1006, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $observed_census_ICU_p_acute_care
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -6.2893, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary
#> 
#> 
#> $N1_N2_avg
#> 
#>  Augmented Dickey-Fuller Test
#> 
#> data:  newX[, i]
#> Dickey-Fuller = -9.3416, Lag order = 7, p-value = 0.01
#> alternative hypothesis: stationary

Further confirmation also comes from the plots of the differenced time series. The trend observed in time series before is now no longer visible. The time series are now looking more stationary.

VARselect has been used from VAR model to determine the most appropriate lag value to regress against each variable in the dataset. VAR model has been run on the differenced time series training data using lag values from 1 to a maximum value of 15 to check the most appropriate lag value. There were 4 error metrics used to determine which lag value from 1 to 15 gave the lowest of those 4 error metrics. In case different error metrics gave different values of lag, then AIC was used as the choice of metric to select the right lag value.

#> $selection
#> AIC(n)  HQ(n)  SC(n) FPE(n) 
#>     15     15     15     15 
#> 
#> $criteria
#>                    1             2             3             4             5
#> AIC(n) -1.216633e+02 -1.219541e+02 -1.230962e+02 -1.235904e+02 -1.266219e+02
#> HQ(n)  -1.199969e+02 -1.186212e+02 -1.180969e+02 -1.169247e+02 -1.182897e+02
#> SC(n)  -1.174663e+02 -1.135600e+02 -1.105050e+02 -1.068021e+02 -1.056366e+02
#> FPE(n)  1.456063e-53  1.104448e-53  3.666841e-54  2.419107e-54  1.332022e-55
#>                    6             7             8             9            10
#> AIC(n) -1.288476e+02 -1.326766e+02 -1.368379e+02 -1.416927e+02 -1.501234e+02
#> HQ(n)  -1.188490e+02 -1.210116e+02 -1.235065e+02 -1.266948e+02 -1.334591e+02
#> SC(n)  -1.036652e+02 -1.032972e+02 -1.032614e+02 -1.039191e+02 -1.081528e+02
#> FPE(n)  1.762924e-56  5.145286e-58  1.210135e-59  1.648414e-61  7.581507e-65
#>                   11            12            13            14            15
#> AIC(n) -1.577191e+02 -1.635819e+02 -1.705205e+02 -1.765209e+02 -1.846040e+02
#> HQ(n)  -1.393884e+02 -1.435848e+02 -1.488570e+02 -1.531909e+02 -1.596076e+02
#> SC(n)  -1.115514e+02 -1.132171e+02 -1.159587e+02 -1.177619e+02 -1.216480e+02
#> FPE(n)  1.023473e-67  1.073714e-69  5.900972e-72  1.518285e-73  1.206351e-75

From the lowest error metric, 15th lag is the most appropriate to use. The VAR model is now fit on the training data with a lag of 15 days.

Summary of the VAR model with beta coefficients for all the 15 different lag values and their statistical significance for all the different variables in the model that are being regressed against each other is available in VAR. Below we can see the summary for the hospitalization variable observed_census_ICU_p_acute_care. We only look at observed_census_ICU_p_acute_care variable since we are only interested in forecasting the hospitalization.

#> 
#> Call:
#> lm(formula = y ~ -1 + ., data = datamat)
#> 
#> Residuals:
#>      Min       1Q   Median       3Q      Max 
#> -0.28900 -0.02439 -0.00056  0.02420  0.33102 
#> 
#> Coefficients:
#>                                             Estimate Std. Error t value
#> Percent_of_Ottawa_residents_5_1_dose.l1    6.516e-02  2.639e-01   0.247
#> Percent_of_Ottawa_residents_12_1_dose.l1  -6.080e-01  2.108e+00  -0.288
#> Percent_of_Ottawa_residents_18_1_dose.l1  -1.176e+00  5.356e+00  -0.220
#> Percent_of_Ottawa_residents_30_1_dose.l1  -5.700e+00  5.480e+00  -1.040
#> Percent_of_Ottawa_residents_40_1_dose.l1   7.750e+00  7.487e+00   1.035
#> Percent_of_Ottawa_residents_50_1_dose.l1  -1.394e+00  6.527e+00  -0.214
#> Percent_of_Ottawa_residents_60_1_dose.l1  -6.012e-01  5.539e+00  -0.109
#> Percent_of_Ottawa_residents_70_1_dose.l1  -2.809e+00  8.481e+00  -0.331
#> Percent_of_Ottawa_residents_80_1_dose.l1  -9.400e+00  1.045e+01  -0.900
#> Percent_of_Ottawa_residents_5_2_dose.l1    3.843e-01  3.110e-01   1.236
#> Percent_of_Ottawa_residents_12_2_dose.l1   3.544e-01  1.020e+00   0.347
#> Percent_of_Ottawa_residents_18_2_dose.l1  -3.749e-02  6.364e-02  -0.589
#> Percent_of_Ottawa_residents_30_2_dose.l1   2.474e+00  4.180e+00   0.592
#> Percent_of_Ottawa_residents_40_2_dose.l1  -3.279e-01  3.239e+00  -0.101
#> Percent_of_Ottawa_residents_50_2_dose.l1  -1.907e+00  3.776e+00  -0.505
#> Percent_of_Ottawa_residents_60_2_dose.l1   8.456e-01  3.887e+00   0.218
#> Percent_of_Ottawa_residents_70_2_dose.l1  -1.211e+00  3.926e+00  -0.309
#> Percent_of_Ottawa_residents_80_2_dose.l1   9.619e+00  4.076e+00   2.360
#> observed_census_ICU_p_acute_care.l1       -9.255e-03  1.357e-01  -0.068
#> N1_N2_avg.l1                              -9.767e-02  2.833e-01  -0.345
#> Percent_of_Ottawa_residents_5_1_dose.l2   -1.284e-01  2.754e-01  -0.466
#> Percent_of_Ottawa_residents_12_1_dose.l2  -6.964e-01  2.051e+00  -0.340
#> Percent_of_Ottawa_residents_18_1_dose.l2  -7.627e-01  4.991e+00  -0.153
#> Percent_of_Ottawa_residents_30_1_dose.l2  -2.012e+00  6.049e+00  -0.333
#> Percent_of_Ottawa_residents_40_1_dose.l2   2.511e+00  6.810e+00   0.369
#> Percent_of_Ottawa_residents_50_1_dose.l2  -8.784e-01  6.370e+00  -0.138
#> Percent_of_Ottawa_residents_60_1_dose.l2  -1.087e+01  5.637e+00  -1.929
#> Percent_of_Ottawa_residents_70_1_dose.l2   6.501e+00  8.815e+00   0.737
#> Percent_of_Ottawa_residents_80_1_dose.l2   1.119e+01  1.152e+01   0.971
#> Percent_of_Ottawa_residents_5_2_dose.l2   -1.827e-02  3.187e-01  -0.057
#> Percent_of_Ottawa_residents_12_2_dose.l2  -4.713e-01  8.093e-01  -0.582
#> Percent_of_Ottawa_residents_18_2_dose.l2   1.607e-02  6.040e-02   0.266
#> Percent_of_Ottawa_residents_30_2_dose.l2   5.430e+00  3.772e+00   1.440
#> Percent_of_Ottawa_residents_40_2_dose.l2   4.085e-01  3.747e+00   0.109
#> Percent_of_Ottawa_residents_50_2_dose.l2   1.107e+00  3.564e+00   0.311
#> Percent_of_Ottawa_residents_60_2_dose.l2   1.784e+00  4.149e+00   0.430
#> Percent_of_Ottawa_residents_70_2_dose.l2  -2.223e+00  4.615e+00  -0.482
#> Percent_of_Ottawa_residents_80_2_dose.l2  -2.251e+00  3.522e+00  -0.639
#> observed_census_ICU_p_acute_care.l2        1.107e-01  1.306e-01   0.847
#> N1_N2_avg.l2                               5.047e-01  3.267e-01   1.545
#> Percent_of_Ottawa_residents_5_1_dose.l3    1.562e-03  2.788e-01   0.006
#> Percent_of_Ottawa_residents_12_1_dose.l3   8.322e-01  1.938e+00   0.429
#> Percent_of_Ottawa_residents_18_1_dose.l3  -6.053e-01  4.831e+00  -0.125
#> Percent_of_Ottawa_residents_30_1_dose.l3   7.471e-01  6.012e+00   0.124
#> Percent_of_Ottawa_residents_40_1_dose.l3  -5.185e+00  6.699e+00  -0.774
#> Percent_of_Ottawa_residents_50_1_dose.l3  -2.697e+00  6.773e+00  -0.398
#> Percent_of_Ottawa_residents_60_1_dose.l3   7.141e-01  5.422e+00   0.132
#> Percent_of_Ottawa_residents_70_1_dose.l3  -6.092e+00  8.445e+00  -0.721
#> Percent_of_Ottawa_residents_80_1_dose.l3   1.247e+00  1.118e+01   0.112
#> Percent_of_Ottawa_residents_5_2_dose.l3    1.238e-02  3.169e-01   0.039
#> Percent_of_Ottawa_residents_12_2_dose.l3   9.575e-01  9.307e-01   1.029
#> Percent_of_Ottawa_residents_18_2_dose.l3   2.908e-02  6.443e-02   0.451
#> Percent_of_Ottawa_residents_30_2_dose.l3   4.168e+00  4.558e+00   0.914
#> Percent_of_Ottawa_residents_40_2_dose.l3  -1.644e+00  3.404e+00  -0.483
#> Percent_of_Ottawa_residents_50_2_dose.l3   3.792e+00  3.062e+00   1.238
#> Percent_of_Ottawa_residents_60_2_dose.l3   3.163e+00  3.319e+00   0.953
#> Percent_of_Ottawa_residents_70_2_dose.l3   2.548e+00  3.907e+00   0.652
#> Percent_of_Ottawa_residents_80_2_dose.l3  -9.354e+00  4.454e+00  -2.100
#> observed_census_ICU_p_acute_care.l3        5.496e-02  1.302e-01   0.422
#> N1_N2_avg.l3                               3.183e-01  3.613e-01   0.881
#> Percent_of_Ottawa_residents_5_1_dose.l4   -2.535e-01  2.753e-01  -0.921
#> Percent_of_Ottawa_residents_12_1_dose.l4   1.037e+00  2.108e+00   0.492
#> Percent_of_Ottawa_residents_18_1_dose.l4   2.452e+00  4.530e+00   0.541
#> Percent_of_Ottawa_residents_30_1_dose.l4   2.574e-02  5.967e+00   0.004
#> Percent_of_Ottawa_residents_40_1_dose.l4  -3.719e+00  6.854e+00  -0.543
#> Percent_of_Ottawa_residents_50_1_dose.l4   3.603e+00  6.735e+00   0.535
#> Percent_of_Ottawa_residents_60_1_dose.l4   1.321e+00  6.145e+00   0.215
#> Percent_of_Ottawa_residents_70_1_dose.l4  -3.819e+00  7.979e+00  -0.479
#> Percent_of_Ottawa_residents_80_1_dose.l4  -9.270e+00  1.067e+01  -0.869
#> Percent_of_Ottawa_residents_5_2_dose.l4   -2.921e-02  3.341e-01  -0.087
#> Percent_of_Ottawa_residents_12_2_dose.l4  -1.344e+00  1.099e+00  -1.224
#> Percent_of_Ottawa_residents_18_2_dose.l4   2.666e-02  6.191e-02   0.431
#> Percent_of_Ottawa_residents_30_2_dose.l4   5.840e+00  4.529e+00   1.289
#> Percent_of_Ottawa_residents_40_2_dose.l4  -1.097e+00  3.726e+00  -0.294
#> Percent_of_Ottawa_residents_50_2_dose.l4   3.198e+00  3.153e+00   1.014
#> Percent_of_Ottawa_residents_60_2_dose.l4   1.298e+00  3.692e+00   0.352
#> Percent_of_Ottawa_residents_70_2_dose.l4   2.185e+00  3.047e+00   0.717
#> Percent_of_Ottawa_residents_80_2_dose.l4  -1.527e+00  5.626e+00  -0.271
#> observed_census_ICU_p_acute_care.l4        8.156e-02  1.272e-01   0.641
#> N1_N2_avg.l4                               3.409e-01  3.569e-01   0.955
#> Percent_of_Ottawa_residents_5_1_dose.l5    1.903e-02  2.868e-01   0.066
#> Percent_of_Ottawa_residents_12_1_dose.l5   3.480e-01  1.926e+00   0.181
#> Percent_of_Ottawa_residents_18_1_dose.l5   5.399e+00  4.579e+00   1.179
#> Percent_of_Ottawa_residents_30_1_dose.l5   9.225e-02  5.698e+00   0.016
#> Percent_of_Ottawa_residents_40_1_dose.l5  -2.864e+00  6.157e+00  -0.465
#> Percent_of_Ottawa_residents_50_1_dose.l5   8.769e+00  6.372e+00   1.376
#> Percent_of_Ottawa_residents_60_1_dose.l5  -6.819e-01  6.493e+00  -0.105
#> Percent_of_Ottawa_residents_70_1_dose.l5   3.311e+00  7.201e+00   0.460
#> Percent_of_Ottawa_residents_80_1_dose.l5  -1.518e+00  8.806e+00  -0.172
#> Percent_of_Ottawa_residents_5_2_dose.l5   -7.742e-02  3.277e-01  -0.236
#> Percent_of_Ottawa_residents_12_2_dose.l5  -2.144e-01  8.876e-01  -0.242
#> Percent_of_Ottawa_residents_18_2_dose.l5  -1.224e-02  6.056e-02  -0.202
#> Percent_of_Ottawa_residents_30_2_dose.l5   3.887e+00  4.580e+00   0.849
#> Percent_of_Ottawa_residents_40_2_dose.l5  -1.499e+00  4.060e+00  -0.369
#> Percent_of_Ottawa_residents_50_2_dose.l5  -3.142e+00  2.937e+00  -1.070
#> Percent_of_Ottawa_residents_60_2_dose.l5  -9.013e-01  4.010e+00  -0.225
#> Percent_of_Ottawa_residents_70_2_dose.l5  -1.352e+00  2.914e+00  -0.464
#> Percent_of_Ottawa_residents_80_2_dose.l5   4.566e+00  5.220e+00   0.875
#> observed_census_ICU_p_acute_care.l5       -1.094e-01  1.213e-01  -0.902
#> N1_N2_avg.l5                               3.222e-02  3.763e-01   0.086
#> Percent_of_Ottawa_residents_5_1_dose.l6   -7.168e-02  2.880e-01  -0.249
#> Percent_of_Ottawa_residents_12_1_dose.l6  -7.104e-01  1.994e+00  -0.356
#> Percent_of_Ottawa_residents_18_1_dose.l6  -5.989e+00  5.640e+00  -1.062
#> Percent_of_Ottawa_residents_30_1_dose.l6   3.572e+00  6.019e+00   0.594
#> Percent_of_Ottawa_residents_40_1_dose.l6   8.307e-01  6.602e+00   0.126
#> Percent_of_Ottawa_residents_50_1_dose.l6   3.781e+00  6.198e+00   0.610
#> Percent_of_Ottawa_residents_60_1_dose.l6  -1.791e+00  6.339e+00  -0.283
#> Percent_of_Ottawa_residents_70_1_dose.l6   4.041e+00  7.628e+00   0.530
#> Percent_of_Ottawa_residents_80_1_dose.l6  -2.148e+00  1.120e+01  -0.192
#> Percent_of_Ottawa_residents_5_2_dose.l6   -1.732e-01  3.214e-01  -0.539
#> Percent_of_Ottawa_residents_12_2_dose.l6  -2.679e-01  8.996e-01  -0.298
#> Percent_of_Ottawa_residents_18_2_dose.l6   1.083e-01  6.267e-02   1.728
#> Percent_of_Ottawa_residents_30_2_dose.l6  -2.375e+00  4.029e+00  -0.589
#> Percent_of_Ottawa_residents_40_2_dose.l6  -2.577e+00  4.359e+00  -0.591
#> Percent_of_Ottawa_residents_50_2_dose.l6  -7.613e+00  3.061e+00  -2.487
#> Percent_of_Ottawa_residents_60_2_dose.l6  -1.727e+00  3.942e+00  -0.438
#> Percent_of_Ottawa_residents_70_2_dose.l6  -1.603e+00  2.928e+00  -0.548
#> Percent_of_Ottawa_residents_80_2_dose.l6   6.183e+00  5.634e+00   1.097
#> observed_census_ICU_p_acute_care.l6       -1.818e-01  1.223e-01  -1.486
#> N1_N2_avg.l6                               1.749e-02  3.900e-01   0.045
#> Percent_of_Ottawa_residents_5_1_dose.l7   -1.560e-01  2.697e-01  -0.578
#> Percent_of_Ottawa_residents_12_1_dose.l7  -2.510e+00  1.955e+00  -1.284
#> Percent_of_Ottawa_residents_18_1_dose.l7  -2.985e+00  5.694e+00  -0.524
#> Percent_of_Ottawa_residents_30_1_dose.l7  -2.456e-01  5.335e+00  -0.046
#> Percent_of_Ottawa_residents_40_1_dose.l7   8.114e+00  6.016e+00   1.349
#> Percent_of_Ottawa_residents_50_1_dose.l7   7.259e+00  6.748e+00   1.076
#> Percent_of_Ottawa_residents_60_1_dose.l7  -2.386e+00  6.396e+00  -0.373
#> Percent_of_Ottawa_residents_70_1_dose.l7  -3.369e-01  6.104e+00  -0.055
#> Percent_of_Ottawa_residents_80_1_dose.l7   3.281e+00  1.067e+01   0.307
#> Percent_of_Ottawa_residents_5_2_dose.l7   -1.096e-01  3.467e-01  -0.316
#> Percent_of_Ottawa_residents_12_2_dose.l7  -6.834e-01  9.221e-01  -0.741
#> Percent_of_Ottawa_residents_18_2_dose.l7  -1.535e-02  6.471e-02  -0.237
#> Percent_of_Ottawa_residents_30_2_dose.l7  -2.386e+00  3.693e+00  -0.646
#> Percent_of_Ottawa_residents_40_2_dose.l7   6.246e-01  4.014e+00   0.156
#> Percent_of_Ottawa_residents_50_2_dose.l7   1.239e+00  3.396e+00   0.365
#> Percent_of_Ottawa_residents_60_2_dose.l7  -1.137e+00  3.613e+00  -0.315
#> Percent_of_Ottawa_residents_70_2_dose.l7   2.568e+00  3.258e+00   0.788
#> Percent_of_Ottawa_residents_80_2_dose.l7  -1.094e+00  4.772e+00  -0.229
#> observed_census_ICU_p_acute_care.l7        1.061e-01  1.282e-01   0.828
#> N1_N2_avg.l7                               1.758e-01  3.875e-01   0.454
#> Percent_of_Ottawa_residents_5_1_dose.l8    1.625e-01  2.651e-01   0.613
#> Percent_of_Ottawa_residents_12_1_dose.l8  -1.461e+00  1.982e+00  -0.737
#> Percent_of_Ottawa_residents_18_1_dose.l8   2.022e+00  5.092e+00   0.397
#> Percent_of_Ottawa_residents_30_1_dose.l8  -6.439e+00  5.726e+00  -1.124
#> Percent_of_Ottawa_residents_40_1_dose.l8   7.220e-01  6.211e+00   0.116
#> Percent_of_Ottawa_residents_50_1_dose.l8  -9.534e-01  6.700e+00  -0.142
#> Percent_of_Ottawa_residents_60_1_dose.l8   4.120e+00  7.447e+00   0.553
#> Percent_of_Ottawa_residents_70_1_dose.l8   1.994e+00  6.491e+00   0.307
#> Percent_of_Ottawa_residents_80_1_dose.l8   3.341e+00  9.927e+00   0.337
#> Percent_of_Ottawa_residents_5_2_dose.l8    2.514e-01  3.666e-01   0.686
#> Percent_of_Ottawa_residents_12_2_dose.l8   2.911e-01  8.536e-01   0.341
#> Percent_of_Ottawa_residents_18_2_dose.l8   3.428e-02  6.896e-02   0.497
#> Percent_of_Ottawa_residents_30_2_dose.l8  -6.116e+00  3.640e+00  -1.680
#> Percent_of_Ottawa_residents_40_2_dose.l8  -2.023e+00  3.359e+00  -0.602
#> Percent_of_Ottawa_residents_50_2_dose.l8   2.001e+00  2.821e+00   0.709
#> Percent_of_Ottawa_residents_60_2_dose.l8  -7.002e+00  4.697e+00  -1.491
#> Percent_of_Ottawa_residents_70_2_dose.l8  -9.334e-02  3.430e+00  -0.027
#> Percent_of_Ottawa_residents_80_2_dose.l8  -1.239e-01  5.388e+00  -0.023
#> observed_census_ICU_p_acute_care.l8        8.301e-02  1.249e-01   0.665
#> N1_N2_avg.l8                               4.145e-01  3.827e-01   1.083
#> Percent_of_Ottawa_residents_5_1_dose.l9   -6.819e-02  2.657e-01  -0.257
#> Percent_of_Ottawa_residents_12_1_dose.l9   1.645e+00  2.035e+00   0.808
#> Percent_of_Ottawa_residents_18_1_dose.l9   5.735e+00  4.798e+00   1.195
#> Percent_of_Ottawa_residents_30_1_dose.l9   3.119e+00  5.519e+00   0.565
#> Percent_of_Ottawa_residents_40_1_dose.l9  -6.314e+00  6.965e+00  -0.907
#> Percent_of_Ottawa_residents_50_1_dose.l9  -4.016e+00  6.802e+00  -0.590
#> Percent_of_Ottawa_residents_60_1_dose.l9  -4.226e-01  7.595e+00  -0.056
#> Percent_of_Ottawa_residents_70_1_dose.l9   1.200e+01  6.568e+00   1.827
#> Percent_of_Ottawa_residents_80_1_dose.l9  -1.088e+00  9.390e+00  -0.116
#> Percent_of_Ottawa_residents_5_2_dose.l9    3.048e-01  3.686e-01   0.827
#> Percent_of_Ottawa_residents_12_2_dose.l9  -5.019e-01  9.713e-01  -0.517
#> Percent_of_Ottawa_residents_18_2_dose.l9  -7.853e-02  6.675e-02  -1.177
#> Percent_of_Ottawa_residents_30_2_dose.l9   8.492e-01  2.976e+00   0.285
#> Percent_of_Ottawa_residents_40_2_dose.l9   2.028e+00  3.423e+00   0.592
#> Percent_of_Ottawa_residents_50_2_dose.l9   5.322e+00  3.325e+00   1.601
#> Percent_of_Ottawa_residents_60_2_dose.l9  -6.073e+00  4.313e+00  -1.408
#> Percent_of_Ottawa_residents_70_2_dose.l9  -3.344e+00  3.584e+00  -0.933
#> Percent_of_Ottawa_residents_80_2_dose.l9  -5.116e-01  4.244e+00  -0.121
#> observed_census_ICU_p_acute_care.l9       -1.071e-01  1.185e-01  -0.904
#> N1_N2_avg.l9                               4.367e-01  3.912e-01   1.116
#> Percent_of_Ottawa_residents_5_1_dose.l10  -1.454e-01  2.645e-01  -0.550
#> Percent_of_Ottawa_residents_12_1_dose.l10 -3.899e-01  1.605e+00  -0.243
#> Percent_of_Ottawa_residents_18_1_dose.l10  3.296e+00  5.197e+00   0.634
#> Percent_of_Ottawa_residents_30_1_dose.l10 -3.805e+00  5.436e+00  -0.700
#> Percent_of_Ottawa_residents_40_1_dose.l10  4.647e+00  6.532e+00   0.711
#> Percent_of_Ottawa_residents_50_1_dose.l10 -2.929e+00  6.694e+00  -0.438
#> Percent_of_Ottawa_residents_60_1_dose.l10  2.067e+00  6.681e+00   0.309
#> Percent_of_Ottawa_residents_70_1_dose.l10 -4.562e-01  5.823e+00  -0.078
#> Percent_of_Ottawa_residents_80_1_dose.l10 -6.819e+00  8.544e+00  -0.798
#> Percent_of_Ottawa_residents_5_2_dose.l10   3.757e-01  3.408e-01   1.102
#> Percent_of_Ottawa_residents_12_2_dose.l10  9.019e-01  8.034e-01   1.123
#> Percent_of_Ottawa_residents_18_2_dose.l10 -7.922e-02  6.539e-02  -1.212
#> Percent_of_Ottawa_residents_30_2_dose.l10  1.969e+00  2.852e+00   0.690
#> Percent_of_Ottawa_residents_40_2_dose.l10  3.031e+00  3.323e+00   0.912
#> Percent_of_Ottawa_residents_50_2_dose.l10  2.501e+00  3.434e+00   0.728
#> Percent_of_Ottawa_residents_60_2_dose.l10 -3.956e+00  4.443e+00  -0.890
#> Percent_of_Ottawa_residents_70_2_dose.l10  3.276e-01  3.482e+00   0.094
#> Percent_of_Ottawa_residents_80_2_dose.l10 -8.058e-01  3.869e+00  -0.208
#> observed_census_ICU_p_acute_care.l10      -1.252e-01  1.218e-01  -1.028
#> N1_N2_avg.l10                              3.630e-01  3.717e-01   0.977
#> Percent_of_Ottawa_residents_5_1_dose.l11   2.386e-01  2.692e-01   0.886
#> Percent_of_Ottawa_residents_12_1_dose.l11 -1.627e+00  1.914e+00  -0.850
#> Percent_of_Ottawa_residents_18_1_dose.l11  2.614e+00  3.781e+00   0.691
#> Percent_of_Ottawa_residents_30_1_dose.l11 -9.523e-01  5.012e+00  -0.190
#> Percent_of_Ottawa_residents_40_1_dose.l11  1.164e+00  6.858e+00   0.170
#> Percent_of_Ottawa_residents_50_1_dose.l11 -4.124e+00  6.355e+00  -0.649
#> Percent_of_Ottawa_residents_60_1_dose.l11  2.915e+00  7.601e+00   0.383
#> Percent_of_Ottawa_residents_70_1_dose.l11 -9.565e+00  6.383e+00  -1.499
#> Percent_of_Ottawa_residents_80_1_dose.l11 -7.342e+00  8.913e+00  -0.824
#> Percent_of_Ottawa_residents_5_2_dose.l11   5.565e-02  3.349e-01   0.166
#> Percent_of_Ottawa_residents_12_2_dose.l11 -1.296e+00  7.462e-01  -1.736
#> Percent_of_Ottawa_residents_18_2_dose.l11  7.788e-04  5.895e-02   0.013
#> Percent_of_Ottawa_residents_30_2_dose.l11  3.442e+00  3.098e+00   1.111
#> Percent_of_Ottawa_residents_40_2_dose.l11  2.019e+00  3.243e+00   0.623
#> Percent_of_Ottawa_residents_50_2_dose.l11 -5.780e+00  3.496e+00  -1.654
#> Percent_of_Ottawa_residents_60_2_dose.l11  3.018e+00  4.518e+00   0.668
#> Percent_of_Ottawa_residents_70_2_dose.l11 -3.112e-02  3.527e+00  -0.009
#> Percent_of_Ottawa_residents_80_2_dose.l11  3.584e+00  4.323e+00   0.829
#> observed_census_ICU_p_acute_care.l11       2.583e-02  1.284e-01   0.201
#> N1_N2_avg.l11                              4.507e-01  3.849e-01   1.171
#> Percent_of_Ottawa_residents_5_1_dose.l12  -1.344e-01  2.685e-01  -0.501
#> Percent_of_Ottawa_residents_12_1_dose.l12  7.793e-01  1.724e+00   0.452
#> Percent_of_Ottawa_residents_18_1_dose.l12 -3.533e+00  4.450e+00  -0.794
#> Percent_of_Ottawa_residents_30_1_dose.l12  5.863e+00  5.018e+00   1.168
#> Percent_of_Ottawa_residents_40_1_dose.l12  1.319e+00  6.417e+00   0.206
#> Percent_of_Ottawa_residents_50_1_dose.l12 -6.769e+00  6.158e+00  -1.099
#> Percent_of_Ottawa_residents_60_1_dose.l12 -8.888e+00  7.288e+00  -1.220
#> Percent_of_Ottawa_residents_70_1_dose.l12 -6.131e+00  7.091e+00  -0.865
#> Percent_of_Ottawa_residents_80_1_dose.l12  7.923e+00  8.183e+00   0.968
#> Percent_of_Ottawa_residents_5_2_dose.l12  -3.096e-01  3.350e-01  -0.924
#> Percent_of_Ottawa_residents_12_2_dose.l12 -6.501e-01  7.743e-01  -0.840
#> Percent_of_Ottawa_residents_18_2_dose.l12  5.737e-03  5.892e-02   0.097
#> Percent_of_Ottawa_residents_30_2_dose.l12  3.620e+00  2.558e+00   1.415
#> Percent_of_Ottawa_residents_40_2_dose.l12  1.535e+00  3.581e+00   0.429
#> Percent_of_Ottawa_residents_50_2_dose.l12 -6.272e+00  3.132e+00  -2.003
#> Percent_of_Ottawa_residents_60_2_dose.l12 -1.837e+00  4.001e+00  -0.459
#> Percent_of_Ottawa_residents_70_2_dose.l12  2.365e+00  3.649e+00   0.648
#> Percent_of_Ottawa_residents_80_2_dose.l12  1.139e+00  3.677e+00   0.310
#> observed_census_ICU_p_acute_care.l12      -1.221e-01  1.224e-01  -0.998
#> N1_N2_avg.l12                              7.620e-02  3.838e-01   0.199
#> Percent_of_Ottawa_residents_5_1_dose.l13   1.245e-01  2.610e-01   0.477
#> Percent_of_Ottawa_residents_12_1_dose.l13 -1.244e-01  1.685e+00  -0.074
#> Percent_of_Ottawa_residents_18_1_dose.l13 -2.195e+00  4.109e+00  -0.534
#> Percent_of_Ottawa_residents_30_1_dose.l13  4.381e+00  4.649e+00   0.942
#> Percent_of_Ottawa_residents_40_1_dose.l13 -5.527e+00  5.634e+00  -0.981
#> Percent_of_Ottawa_residents_50_1_dose.l13  5.907e-01  6.337e+00   0.093
#> Percent_of_Ottawa_residents_60_1_dose.l13  6.031e+00  7.083e+00   0.852
#> Percent_of_Ottawa_residents_70_1_dose.l13  5.515e+00  6.100e+00   0.904
#> Percent_of_Ottawa_residents_80_1_dose.l13  9.262e+00  8.205e+00   1.129
#> Percent_of_Ottawa_residents_5_2_dose.l13   1.758e-02  3.248e-01   0.054
#> Percent_of_Ottawa_residents_12_2_dose.l13 -3.745e-01  7.282e-01  -0.514
#> Percent_of_Ottawa_residents_18_2_dose.l13 -3.821e-02  5.667e-02  -0.674
#> Percent_of_Ottawa_residents_30_2_dose.l13  1.102e+00  2.508e+00   0.440
#> Percent_of_Ottawa_residents_40_2_dose.l13  9.512e-02  3.371e+00   0.028
#> Percent_of_Ottawa_residents_50_2_dose.l13 -6.125e+00  3.273e+00  -1.871
#> Percent_of_Ottawa_residents_60_2_dose.l13 -2.534e+00  3.482e+00  -0.728
#> Percent_of_Ottawa_residents_70_2_dose.l13  2.168e+00  4.025e+00   0.539
#> Percent_of_Ottawa_residents_80_2_dose.l13 -1.488e+00  3.913e+00  -0.380
#> observed_census_ICU_p_acute_care.l13      -2.738e-02  1.221e-01  -0.224
#> N1_N2_avg.l13                              7.833e-02  3.982e-01   0.197
#> Percent_of_Ottawa_residents_5_1_dose.l14  -1.045e-03  2.498e-01  -0.004
#> Percent_of_Ottawa_residents_12_1_dose.l14  1.400e+00  1.899e+00   0.737
#> Percent_of_Ottawa_residents_18_1_dose.l14 -2.140e+00  3.838e+00  -0.558
#> Percent_of_Ottawa_residents_30_1_dose.l14  1.963e+00  4.139e+00   0.474
#> Percent_of_Ottawa_residents_40_1_dose.l14 -8.445e+00  5.481e+00  -1.541
#> Percent_of_Ottawa_residents_50_1_dose.l14 -4.561e+00  6.428e+00  -0.710
#> Percent_of_Ottawa_residents_60_1_dose.l14  8.064e+00  7.414e+00   1.088
#> Percent_of_Ottawa_residents_70_1_dose.l14 -3.301e-01  5.252e+00  -0.063
#> Percent_of_Ottawa_residents_80_1_dose.l14  2.080e+00  6.484e+00   0.321
#> Percent_of_Ottawa_residents_5_2_dose.l14  -1.117e-02  3.689e-01  -0.030
#> Percent_of_Ottawa_residents_12_2_dose.l14 -1.297e-01  7.871e-01  -0.165
#> Percent_of_Ottawa_residents_18_2_dose.l14  1.142e-01  5.171e-02   2.209
#> Percent_of_Ottawa_residents_30_2_dose.l14 -2.113e+00  2.492e+00  -0.848
#> Percent_of_Ottawa_residents_40_2_dose.l14 -8.232e-01  3.680e+00  -0.224
#> Percent_of_Ottawa_residents_50_2_dose.l14  3.113e+00  4.328e+00   0.719
#> Percent_of_Ottawa_residents_60_2_dose.l14 -6.688e-01  3.366e+00  -0.199
#> Percent_of_Ottawa_residents_70_2_dose.l14  1.179e+00  3.615e+00   0.326
#> Percent_of_Ottawa_residents_80_2_dose.l14  7.037e-01  3.538e+00   0.199
#> observed_census_ICU_p_acute_care.l14      -3.141e-02  1.144e-01  -0.274
#> N1_N2_avg.l14                              2.169e-01  3.848e-01   0.564
#> Percent_of_Ottawa_residents_5_1_dose.l15  -1.967e-01  2.488e-01  -0.791
#> Percent_of_Ottawa_residents_12_1_dose.l15  3.685e-01  1.851e+00   0.199
#> Percent_of_Ottawa_residents_18_1_dose.l15 -1.034e+00  3.761e+00  -0.275
#> Percent_of_Ottawa_residents_30_1_dose.l15  3.766e+00  4.015e+00   0.938
#> Percent_of_Ottawa_residents_40_1_dose.l15 -1.610e+00  5.368e+00  -0.300
#> Percent_of_Ottawa_residents_50_1_dose.l15 -1.884e+00  5.415e+00  -0.348
#> Percent_of_Ottawa_residents_60_1_dose.l15  3.628e+00  6.632e+00   0.547
#> Percent_of_Ottawa_residents_70_1_dose.l15  1.684e+00  4.744e+00   0.355
#> Percent_of_Ottawa_residents_80_1_dose.l15 -9.497e+00  7.123e+00  -1.333
#> Percent_of_Ottawa_residents_5_2_dose.l15  -4.075e-01  3.819e-01  -1.067
#> Percent_of_Ottawa_residents_12_2_dose.l15  6.979e-01  7.409e-01   0.942
#> Percent_of_Ottawa_residents_18_2_dose.l15  1.258e-02  5.321e-02   0.236
#> Percent_of_Ottawa_residents_30_2_dose.l15 -7.593e-01  2.285e+00  -0.332
#> Percent_of_Ottawa_residents_40_2_dose.l15 -3.210e+00  3.327e+00  -0.965
#> Percent_of_Ottawa_residents_50_2_dose.l15  6.371e+00  3.803e+00   1.675
#> Percent_of_Ottawa_residents_60_2_dose.l15 -3.460e+00  3.133e+00  -1.104
#> Percent_of_Ottawa_residents_70_2_dose.l15  1.056e+00  3.326e+00   0.317
#> Percent_of_Ottawa_residents_80_2_dose.l15  4.499e-01  3.038e+00   0.148
#> observed_census_ICU_p_acute_care.l15       6.331e-02  1.162e-01   0.545
#> N1_N2_avg.l15                             -1.333e-01  3.173e-01  -0.420
#>                                           Pr(>|t|)  
#> Percent_of_Ottawa_residents_5_1_dose.l1     0.8056  
#> Percent_of_Ottawa_residents_12_1_dose.l1    0.7738  
#> Percent_of_Ottawa_residents_18_1_dose.l1    0.8268  
#> Percent_of_Ottawa_residents_30_1_dose.l1    0.3016  
#> Percent_of_Ottawa_residents_40_1_dose.l1    0.3040  
#> Percent_of_Ottawa_residents_50_1_dose.l1    0.8315  
#> Percent_of_Ottawa_residents_60_1_dose.l1    0.9139  
#> Percent_of_Ottawa_residents_70_1_dose.l1    0.7415  
#> Percent_of_Ottawa_residents_80_1_dose.l1    0.3712  
#> Percent_of_Ottawa_residents_5_2_dose.l1     0.2204  
#> Percent_of_Ottawa_residents_12_2_dose.l1    0.7292  
#> Percent_of_Ottawa_residents_18_2_dose.l1    0.5576  
#> Percent_of_Ottawa_residents_30_2_dose.l1    0.5558  
#> Percent_of_Ottawa_residents_40_2_dose.l1    0.9196  
#> Percent_of_Ottawa_residents_50_2_dose.l1    0.6150  
#> Percent_of_Ottawa_residents_60_2_dose.l1    0.8284  
#> Percent_of_Ottawa_residents_70_2_dose.l1    0.7586  
#> Percent_of_Ottawa_residents_80_2_dose.l1    0.0209 *
#> observed_census_ICU_p_acute_care.l1         0.9458  
#> N1_N2_avg.l1                                0.7313  
#> Percent_of_Ottawa_residents_5_1_dose.l2     0.6424  
#> Percent_of_Ottawa_residents_12_1_dose.l2    0.7351  
#> Percent_of_Ottawa_residents_18_1_dose.l2    0.8789  
#> Percent_of_Ottawa_residents_30_1_dose.l2    0.7404  
#> Percent_of_Ottawa_residents_40_1_dose.l2    0.7134  
#> Percent_of_Ottawa_residents_50_1_dose.l2    0.8907  
#> Percent_of_Ottawa_residents_60_1_dose.l2    0.0576 .
#> Percent_of_Ottawa_residents_70_1_dose.l2    0.4632  
#> Percent_of_Ottawa_residents_80_1_dose.l2    0.3345  
#> Percent_of_Ottawa_residents_5_2_dose.l2     0.9544  
#> Percent_of_Ottawa_residents_12_2_dose.l2    0.5621  
#> Percent_of_Ottawa_residents_18_2_dose.l2    0.7910  
#> Percent_of_Ottawa_residents_30_2_dose.l2    0.1542  
#> Percent_of_Ottawa_residents_40_2_dose.l2    0.9135  
#> Percent_of_Ottawa_residents_50_2_dose.l2    0.7570  
#> Percent_of_Ottawa_residents_60_2_dose.l2    0.6684  
#> Percent_of_Ottawa_residents_70_2_dose.l2    0.6315  
#> Percent_of_Ottawa_residents_80_2_dose.l2    0.5247  
#> observed_census_ICU_p_acute_care.l2         0.3995  
#> N1_N2_avg.l2                                0.1267  
#> Percent_of_Ottawa_residents_5_1_dose.l3     0.9955  
#> Percent_of_Ottawa_residents_12_1_dose.l3    0.6688  
#> Percent_of_Ottawa_residents_18_1_dose.l3    0.9006  
#> Percent_of_Ottawa_residents_30_1_dose.l3    0.9014  
#> Percent_of_Ottawa_residents_40_1_dose.l3    0.4414  
#> Percent_of_Ottawa_residents_50_1_dose.l3    0.6916  
#> Percent_of_Ottawa_residents_60_1_dose.l3    0.8956  
#> Percent_of_Ottawa_residents_70_1_dose.l3    0.4730  
#> Percent_of_Ottawa_residents_80_1_dose.l3    0.9115  
#> Percent_of_Ottawa_residents_5_2_dose.l3     0.9689  
#> Percent_of_Ottawa_residents_12_2_dose.l3    0.3069  
#> Percent_of_Ottawa_residents_18_2_dose.l3    0.6531  
#> Percent_of_Ottawa_residents_30_2_dose.l3    0.3635  
#> Percent_of_Ottawa_residents_40_2_dose.l3    0.6306  
#> Percent_of_Ottawa_residents_50_2_dose.l3    0.2195  
#> Percent_of_Ottawa_residents_60_2_dose.l3    0.3437  
#> Percent_of_Ottawa_residents_70_2_dose.l3    0.5163  
#> Percent_of_Ottawa_residents_80_2_dose.l3    0.0391 *
#> observed_census_ICU_p_acute_care.l3         0.6742  
#> N1_N2_avg.l3                                0.3810  
#> Percent_of_Ottawa_residents_5_1_dose.l4     0.3601  
#> Percent_of_Ottawa_residents_12_1_dose.l4    0.6241  
#> Percent_of_Ottawa_residents_18_1_dose.l4    0.5899  
#> Percent_of_Ottawa_residents_30_1_dose.l4    0.9966  
#> Percent_of_Ottawa_residents_40_1_dose.l4    0.5891  
#> Percent_of_Ottawa_residents_50_1_dose.l4    0.5943  
#> Percent_of_Ottawa_residents_60_1_dose.l4    0.8304  
#> Percent_of_Ottawa_residents_70_1_dose.l4    0.6336  
#> Percent_of_Ottawa_residents_80_1_dose.l4    0.3878  
#> Percent_of_Ottawa_residents_5_2_dose.l4     0.9306  
#> Percent_of_Ottawa_residents_12_2_dose.l4    0.2250  
#> Percent_of_Ottawa_residents_18_2_dose.l4    0.6680  
#> Percent_of_Ottawa_residents_30_2_dose.l4    0.2013  
#> Percent_of_Ottawa_residents_40_2_dose.l4    0.7693  
#> Percent_of_Ottawa_residents_50_2_dose.l4    0.3137  
#> Percent_of_Ottawa_residents_60_2_dose.l4    0.7260  
#> Percent_of_Ottawa_residents_70_2_dose.l4    0.4756  
#> Percent_of_Ottawa_residents_80_2_dose.l4    0.7869  
#> observed_census_ICU_p_acute_care.l4         0.5232  
#> N1_N2_avg.l4                                0.3425  
#> Percent_of_Ottawa_residents_5_1_dose.l5     0.9473  
#> Percent_of_Ottawa_residents_12_1_dose.l5    0.8571  
#> Percent_of_Ottawa_residents_18_1_dose.l5    0.2421  
#> Percent_of_Ottawa_residents_30_1_dose.l5    0.9871  
#> Percent_of_Ottawa_residents_40_1_dose.l5    0.6432  
#> Percent_of_Ottawa_residents_50_1_dose.l5    0.1729  
#> Percent_of_Ottawa_residents_60_1_dose.l5    0.9166  
#> Percent_of_Ottawa_residents_70_1_dose.l5    0.6470  
#> Percent_of_Ottawa_residents_80_1_dose.l5    0.8636  
#> Percent_of_Ottawa_residents_5_2_dose.l5     0.8139  
#> Percent_of_Ottawa_residents_12_2_dose.l5    0.8098  
#> Percent_of_Ottawa_residents_18_2_dose.l5    0.8403  
#> Percent_of_Ottawa_residents_30_2_dose.l5    0.3989  
#> Percent_of_Ottawa_residents_40_2_dose.l5    0.7131  
#> Percent_of_Ottawa_residents_50_2_dose.l5    0.2881  
#> Percent_of_Ottawa_residents_60_2_dose.l5    0.8228  
#> Percent_of_Ottawa_residents_70_2_dose.l5    0.6442  
#> Percent_of_Ottawa_residents_80_2_dose.l5    0.3845  
#> observed_census_ICU_p_acute_care.l5         0.3701  
#> N1_N2_avg.l5                                0.9320  
#> Percent_of_Ottawa_residents_5_1_dose.l6     0.8041  
#> Percent_of_Ottawa_residents_12_1_dose.l6    0.7227  
#> Percent_of_Ottawa_residents_18_1_dose.l6    0.2917  
#> Percent_of_Ottawa_residents_30_1_dose.l6    0.5546  
#> Percent_of_Ottawa_residents_40_1_dose.l6    0.9002  
#> Percent_of_Ottawa_residents_50_1_dose.l6    0.5436  
#> Percent_of_Ottawa_residents_60_1_dose.l6    0.7783  
#> Percent_of_Ottawa_residents_70_1_dose.l6    0.5978  
#> Percent_of_Ottawa_residents_80_1_dose.l6    0.8484  
#> Percent_of_Ottawa_residents_5_2_dose.l6     0.5916  
#> Percent_of_Ottawa_residents_12_2_dose.l6    0.7667  
#> Percent_of_Ottawa_residents_18_2_dose.l6    0.0882 .
#> Percent_of_Ottawa_residents_30_2_dose.l6    0.5574  
#> Percent_of_Ottawa_residents_40_2_dose.l6    0.5562  
#> Percent_of_Ottawa_residents_50_2_dose.l6    0.0152 *
#> Percent_of_Ottawa_residents_60_2_dose.l6    0.6626  
#> Percent_of_Ottawa_residents_70_2_dose.l6    0.5857  
#> Percent_of_Ottawa_residents_80_2_dose.l6    0.2760  
#> observed_census_ICU_p_acute_care.l6         0.1415  
#> N1_N2_avg.l6                                0.9643  
#> Percent_of_Ottawa_residents_5_1_dose.l7     0.5649  
#> Percent_of_Ottawa_residents_12_1_dose.l7    0.2031  
#> Percent_of_Ottawa_residents_18_1_dose.l7    0.6017  
#> Percent_of_Ottawa_residents_30_1_dose.l7    0.9634  
#> Percent_of_Ottawa_residents_40_1_dose.l7    0.1816  
#> Percent_of_Ottawa_residents_50_1_dose.l7    0.2856  
#> Percent_of_Ottawa_residents_60_1_dose.l7    0.7102  
#> Percent_of_Ottawa_residents_70_1_dose.l7    0.9561  
#> Percent_of_Ottawa_residents_80_1_dose.l7    0.7594  
#> Percent_of_Ottawa_residents_5_2_dose.l7     0.7529  
#> Percent_of_Ottawa_residents_12_2_dose.l7    0.4610  
#> Percent_of_Ottawa_residents_18_2_dose.l7    0.8132  
#> Percent_of_Ottawa_residents_30_2_dose.l7    0.5203  
#> Percent_of_Ottawa_residents_40_2_dose.l7    0.8768  
#> Percent_of_Ottawa_residents_50_2_dose.l7    0.7162  
#> Percent_of_Ottawa_residents_60_2_dose.l7    0.7539  
#> Percent_of_Ottawa_residents_70_2_dose.l7    0.4330  
#> Percent_of_Ottawa_residents_80_2_dose.l7    0.8193  
#> observed_census_ICU_p_acute_care.l7         0.4104  
#> N1_N2_avg.l7                                0.6514  
#> Percent_of_Ottawa_residents_5_1_dose.l8     0.5418  
#> Percent_of_Ottawa_residents_12_1_dose.l8    0.4633  
#> Percent_of_Ottawa_residents_18_1_dose.l8    0.6924  
#> Percent_of_Ottawa_residents_30_1_dose.l8    0.2645  
#> Percent_of_Ottawa_residents_40_1_dose.l8    0.9078  
#> Percent_of_Ottawa_residents_50_1_dose.l8    0.8872  
#> Percent_of_Ottawa_residents_60_1_dose.l8    0.5818  
#> Percent_of_Ottawa_residents_70_1_dose.l8    0.7596  
#> Percent_of_Ottawa_residents_80_1_dose.l8    0.7374  
#> Percent_of_Ottawa_residents_5_2_dose.l8     0.4950  
#> Percent_of_Ottawa_residents_12_2_dose.l8    0.7341  
#> Percent_of_Ottawa_residents_18_2_dose.l8    0.6206  
#> Percent_of_Ottawa_residents_30_2_dose.l8    0.0971 .
#> Percent_of_Ottawa_residents_40_2_dose.l8    0.5488  
#> Percent_of_Ottawa_residents_50_2_dose.l8    0.4802  
#> Percent_of_Ottawa_residents_60_2_dose.l8    0.1403  
#> Percent_of_Ottawa_residents_70_2_dose.l8    0.9784  
#> Percent_of_Ottawa_residents_80_2_dose.l8    0.9817  
#> observed_census_ICU_p_acute_care.l8         0.5084  
#> N1_N2_avg.l8                                0.2822  
#> Percent_of_Ottawa_residents_5_1_dose.l9     0.7982  
#> Percent_of_Ottawa_residents_12_1_dose.l9    0.4215  
#> Percent_of_Ottawa_residents_18_1_dose.l9    0.2357  
#> Percent_of_Ottawa_residents_30_1_dose.l9    0.5738  
#> Percent_of_Ottawa_residents_40_1_dose.l9    0.3676  
#> Percent_of_Ottawa_residents_50_1_dose.l9    0.5567  
#> Percent_of_Ottawa_residents_60_1_dose.l9    0.9558  
#> Percent_of_Ottawa_residents_70_1_dose.l9    0.0718 .
#> Percent_of_Ottawa_residents_80_1_dose.l9    0.9081  
#> Percent_of_Ottawa_residents_5_2_dose.l9     0.4109  
#> Percent_of_Ottawa_residents_12_2_dose.l9    0.6068  
#> Percent_of_Ottawa_residents_18_2_dose.l9    0.2431  
#> Percent_of_Ottawa_residents_30_2_dose.l9    0.7762  
#> Percent_of_Ottawa_residents_40_2_dose.l9    0.5554  
#> Percent_of_Ottawa_residents_50_2_dose.l9    0.1137  
#> Percent_of_Ottawa_residents_60_2_dose.l9    0.1633  
#> Percent_of_Ottawa_residents_70_2_dose.l9    0.3538  
#> Percent_of_Ottawa_residents_80_2_dose.l9    0.9044  
#> observed_census_ICU_p_acute_care.l9         0.3687  
#> N1_N2_avg.l9                                0.2679  
#> Percent_of_Ottawa_residents_5_1_dose.l10    0.5840  
#> Percent_of_Ottawa_residents_12_1_dose.l10   0.8087  
#> Percent_of_Ottawa_residents_18_1_dose.l10   0.5279  
#> Percent_of_Ottawa_residents_30_1_dose.l10   0.4861  
#> Percent_of_Ottawa_residents_40_1_dose.l10   0.4790  
#> Percent_of_Ottawa_residents_50_1_dose.l10   0.6630  
#> Percent_of_Ottawa_residents_60_1_dose.l10   0.7580  
#> Percent_of_Ottawa_residents_70_1_dose.l10   0.9378  
#> Percent_of_Ottawa_residents_80_1_dose.l10   0.4274  
#> Percent_of_Ottawa_residents_5_2_dose.l10    0.2739  
#> Percent_of_Ottawa_residents_12_2_dose.l10   0.2652  
#> Percent_of_Ottawa_residents_18_2_dose.l10   0.2295  
#> Percent_of_Ottawa_residents_30_2_dose.l10   0.4921  
#> Percent_of_Ottawa_residents_40_2_dose.l10   0.3647  
#> Percent_of_Ottawa_residents_50_2_dose.l10   0.4686  
#> Percent_of_Ottawa_residents_60_2_dose.l10   0.3761  
#> Percent_of_Ottawa_residents_70_2_dose.l10   0.9253  
#> Percent_of_Ottawa_residents_80_2_dose.l10   0.8356  
#> observed_census_ICU_p_acute_care.l10        0.3071  
#> N1_N2_avg.l10                               0.3319  
#> Percent_of_Ottawa_residents_5_1_dose.l11    0.3783  
#> Percent_of_Ottawa_residents_12_1_dose.l11   0.3980  
#> Percent_of_Ottawa_residents_18_1_dose.l11   0.4914  
#> Percent_of_Ottawa_residents_30_1_dose.l11   0.8498  
#> Percent_of_Ottawa_residents_40_1_dose.l11   0.8657  
#> Percent_of_Ottawa_residents_50_1_dose.l11   0.5184  
#> Percent_of_Ottawa_residents_60_1_dose.l11   0.7025  
#> Percent_of_Ottawa_residents_70_1_dose.l11   0.1382  
#> Percent_of_Ottawa_residents_80_1_dose.l11   0.4127  
#> Percent_of_Ottawa_residents_5_2_dose.l11    0.8685  
#> Percent_of_Ottawa_residents_12_2_dose.l11   0.0866 .
#> Percent_of_Ottawa_residents_18_2_dose.l11   0.9895  
#> Percent_of_Ottawa_residents_30_2_dose.l11   0.2701  
#> Percent_of_Ottawa_residents_40_2_dose.l11   0.5355  
#> Percent_of_Ottawa_residents_50_2_dose.l11   0.1025  
#> Percent_of_Ottawa_residents_60_2_dose.l11   0.5061  
#> Percent_of_Ottawa_residents_70_2_dose.l11   0.9930  
#> Percent_of_Ottawa_residents_80_2_dose.l11   0.4097  
#> observed_census_ICU_p_acute_care.l11        0.8411  
#> N1_N2_avg.l11                               0.2453  
#> Percent_of_Ottawa_residents_5_1_dose.l12    0.6181  
#> Percent_of_Ottawa_residents_12_1_dose.l12   0.6525  
#> Percent_of_Ottawa_residents_18_1_dose.l12   0.4297  
#> Percent_of_Ottawa_residents_30_1_dose.l12   0.2464  
#> Percent_of_Ottawa_residents_40_1_dose.l12   0.8377  
#> Percent_of_Ottawa_residents_50_1_dose.l12   0.2753  
#> Percent_of_Ottawa_residents_60_1_dose.l12   0.2265  
#> Percent_of_Ottawa_residents_70_1_dose.l12   0.3901  
#> Percent_of_Ottawa_residents_80_1_dose.l12   0.3361  
#> Percent_of_Ottawa_residents_5_2_dose.l12    0.3583  
#> Percent_of_Ottawa_residents_12_2_dose.l12   0.4038  
#> Percent_of_Ottawa_residents_18_2_dose.l12   0.9227  
#> Percent_of_Ottawa_residents_30_2_dose.l12   0.1613  
#> Percent_of_Ottawa_residents_40_2_dose.l12   0.6694  
#> Percent_of_Ottawa_residents_50_2_dose.l12   0.0489 *
#> Percent_of_Ottawa_residents_60_2_dose.l12   0.6475  
#> Percent_of_Ottawa_residents_70_2_dose.l12   0.5188  
#> Percent_of_Ottawa_residents_80_2_dose.l12   0.7576  
#> observed_census_ICU_p_acute_care.l12        0.3216  
#> N1_N2_avg.l12                               0.8432  
#> Percent_of_Ottawa_residents_5_1_dose.l13    0.6346  
#> Percent_of_Ottawa_residents_12_1_dose.l13   0.9414  
#> Percent_of_Ottawa_residents_18_1_dose.l13   0.5948  
#> Percent_of_Ottawa_residents_30_1_dose.l13   0.3491  
#> Percent_of_Ottawa_residents_40_1_dose.l13   0.3298  
#> Percent_of_Ottawa_residents_50_1_dose.l13   0.9260  
#> Percent_of_Ottawa_residents_60_1_dose.l13   0.3972  
#> Percent_of_Ottawa_residents_70_1_dose.l13   0.3689  
#> Percent_of_Ottawa_residents_80_1_dose.l13   0.2626  
#> Percent_of_Ottawa_residents_5_2_dose.l13    0.9570  
#> Percent_of_Ottawa_residents_12_2_dose.l13   0.6086  
#> Percent_of_Ottawa_residents_18_2_dose.l13   0.5023  
#> Percent_of_Ottawa_residents_30_2_dose.l13   0.6616  
#> Percent_of_Ottawa_residents_40_2_dose.l13   0.9776  
#> Percent_of_Ottawa_residents_50_2_dose.l13   0.0653 .
#> Percent_of_Ottawa_residents_60_2_dose.l13   0.4691  
#> Percent_of_Ottawa_residents_70_2_dose.l13   0.5917  
#> Percent_of_Ottawa_residents_80_2_dose.l13   0.7049  
#> observed_census_ICU_p_acute_care.l13        0.8232  
#> N1_N2_avg.l13                               0.8446  
#> Percent_of_Ottawa_residents_5_1_dose.l14    0.9967  
#> Percent_of_Ottawa_residents_12_1_dose.l14   0.4633  
#> Percent_of_Ottawa_residents_18_1_dose.l14   0.5788  
#> Percent_of_Ottawa_residents_30_1_dose.l14   0.6367  
#> Percent_of_Ottawa_residents_40_1_dose.l14   0.1276  
#> Percent_of_Ottawa_residents_50_1_dose.l14   0.4802  
#> Percent_of_Ottawa_residents_60_1_dose.l14   0.2803  
#> Percent_of_Ottawa_residents_70_1_dose.l14   0.9500  
#> Percent_of_Ottawa_residents_80_1_dose.l14   0.7493  
#> Percent_of_Ottawa_residents_5_2_dose.l14    0.9759  
#> Percent_of_Ottawa_residents_12_2_dose.l14   0.8695  
#> Percent_of_Ottawa_residents_18_2_dose.l14   0.0303 *
#> Percent_of_Ottawa_residents_30_2_dose.l14   0.3991  
#> Percent_of_Ottawa_residents_40_2_dose.l14   0.8236  
#> Percent_of_Ottawa_residents_50_2_dose.l14   0.4743  
#> Percent_of_Ottawa_residents_60_2_dose.l14   0.8431  
#> Percent_of_Ottawa_residents_70_2_dose.l14   0.7452  
#> Percent_of_Ottawa_residents_80_2_dose.l14   0.8429  
#> observed_census_ICU_p_acute_care.l14        0.7845  
#> N1_N2_avg.l14                               0.5746  
#> Percent_of_Ottawa_residents_5_1_dose.l15    0.4316  
#> Percent_of_Ottawa_residents_12_1_dose.l15   0.8427  
#> Percent_of_Ottawa_residents_18_1_dose.l15   0.7841  
#> Percent_of_Ottawa_residents_30_1_dose.l15   0.3513  
#> Percent_of_Ottawa_residents_40_1_dose.l15   0.7650  
#> Percent_of_Ottawa_residents_50_1_dose.l15   0.7289  
#> Percent_of_Ottawa_residents_60_1_dose.l15   0.5860  
#> Percent_of_Ottawa_residents_70_1_dose.l15   0.7237  
#> Percent_of_Ottawa_residents_80_1_dose.l15   0.1865  
#> Percent_of_Ottawa_residents_5_2_dose.l15    0.2894  
#> Percent_of_Ottawa_residents_12_2_dose.l15   0.3493  
#> Percent_of_Ottawa_residents_18_2_dose.l15   0.8137  
#> Percent_of_Ottawa_residents_30_2_dose.l15   0.7406  
#> Percent_of_Ottawa_residents_40_2_dose.l15   0.3378  
#> Percent_of_Ottawa_residents_50_2_dose.l15   0.0981 .
#> Percent_of_Ottawa_residents_60_2_dose.l15   0.2730  
#> Percent_of_Ottawa_residents_70_2_dose.l15   0.7518  
#> Percent_of_Ottawa_residents_80_2_dose.l15   0.8827  
#> observed_census_ICU_p_acute_care.l15        0.5875  
#> N1_N2_avg.l15                               0.6756  
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#> 
#> Residual standard error: 0.1963 on 74 degrees of freedom
#> Multiple R-squared:  0.7117, Adjusted R-squared:  -0.4572 
#> F-statistic: 0.6089 on 300 and 74 DF,  p-value: 0.9979

SERIAL CORRELATION BETWEEN ERRORS

The serial correlation between errors have been tested using Portmanteau Test (asymototic). The null hypothesis of this test states that there is no correlation between errors. The p-value being statistically significant shows that there is correlation between the errors. The model may need further tweaking for reducing the correlation between errors.

#> 
#>  Portmanteau Test (asymptotic)
#> 
#> data:  Residuals of VAR object var.a
#> Chi-squared = 12060, df = 400, p-value < 2.2e-16

Predicting 6 days ahead.

FORECASTING 6 DAYS AHEAD INTO FUTURE and CALCULATING ERROR ON THE TEST DATA SET

The difference has been readded back to the predictions to bring it back on original scale.

RMSE Error metric on test data

#> [1] 6.500961

MAPE Error metric on test data

#> [1] 0.06514741

Plots comparing actual test data and predicted test data

Summary:

In Summary, the error metrics on the test data looks reasonably good with a low RMSE of 6.5 compared to the standard deviation of original data which is about 35. The MAPE metric is also around 6% which is quite good. For predictive analysis VAR model can be very effective. Further tuning of hyperparameters can be experimented on this model for better performance.

VAR model is based on the fact that one timeseries causes changes in other time series in multivariate data.

The Granger causality test has been performed to confirm if a timeseries from our data causes another timeseries which means if the lagged values of the given time series being tested is causing the other time series in the data.

The Null Hypothesis H0 is: The given time series does not cause the other time series in the data. Alternate hypthesis H1: The given time series causes the other time series in the data.

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_5_1_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_5_1_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 0.64794, df1 = 285, df2 = 1480, p-value = 1
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_5_1_dose and
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 49.842, df = 19, p-value = 0.0001383

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_5_2_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_5_2_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 1.2235, df1 = 285, df2 = 1480, p-value = 0.01146
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_5_2_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 86.532, df = 19, p-value = 1.357e-10

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_12_1_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_12_1_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.5232, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_12_1_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 106.57, df = 19, p-value = 3.397e-14

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_12_2_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_12_2_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.3088, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_12_2_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 116.1, df = 19, p-value = 5.551e-16

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_18_1_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_18_1_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.0602, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_18_1_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 122.87, df = 19, p-value < 2.2e-16

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_18_2_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_18_2_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 1.4371, df1 = 285, df2 = 1480, p-value = 1.697e-05
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_18_2_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 82.769, df = 19, p-value = 6.167e-10

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_30_1_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_30_1_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 2.647, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_30_1_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 70.221, df = 19, p-value = 8.451e-08

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_30_2_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_30_2_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.8723, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_30_2_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 104.81, df = 19, p-value = 7.139e-14

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_40_1_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_40_1_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.2973, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_40_1_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 124.64, df = 19, p-value < 2.2e-16

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_40_2_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_40_2_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.3536, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_40_2_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 128.73, df = 19, p-value < 2.2e-16

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_50_1_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_50_1_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.1835, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_50_1_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 87.226, df = 19, p-value = 1.025e-10

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_50_2_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_50_2_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 4.3296, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_50_2_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 129.89, df = 19, p-value < 2.2e-16

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_60_1_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_60_1_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.2732, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_60_1_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 93.7, df = 19, p-value = 7.284e-12

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_60_2_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_60_2_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.1245, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_60_2_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 113.32, df = 19, p-value = 1.998e-15

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_70_1_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_70_1_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 4.1005, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_70_1_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 127.23, df = 19, p-value < 2.2e-16

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_70_2_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_70_2_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 3.251, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_70_2_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 61.962, df = 19, p-value = 1.886e-06

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_80_1_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_80_1_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 4.2452, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_80_1_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 124.51, df = 19, p-value < 2.2e-16

From the statistically significant p-value, it shows that Percent_of_Ottawa_residents_80_2_dose causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: Percent_of_Ottawa_residents_80_2_dose do not
#>  Granger-cause Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> F-Test = 4.0135, df1 = 285, df2 = 1480, p-value < 2.2e-16
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between:
#>  Percent_of_Ottawa_residents_80_2_dose and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose observed_census_ICU_p_acute_care
#>  N1_N2_avg
#> 
#> data:  VAR object var.a
#> Chi-squared = 91.573, df = 19, p-value = 1.744e-11

From the statistically significant p-value, it shows that wastewater signal N1_N2_avg causes the other time series in the data.

#> $Granger
#> 
#>  Granger causality H0: N1_N2_avg do not Granger-cause
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#> 
#> data:  VAR object var.a
#> F-Test = 1.1329, df1 = 285, df2 = 1480, p-value = 0.0803
#> 
#> 
#> $Instant
#> 
#>  H0: No instantaneous causality between: N1_N2_avg and
#>  Percent_of_Ottawa_residents_5_1_dose
#>  Percent_of_Ottawa_residents_12_1_dose
#>  Percent_of_Ottawa_residents_18_1_dose
#>  Percent_of_Ottawa_residents_30_1_dose
#>  Percent_of_Ottawa_residents_40_1_dose
#>  Percent_of_Ottawa_residents_50_1_dose
#>  Percent_of_Ottawa_residents_60_1_dose
#>  Percent_of_Ottawa_residents_70_1_dose
#>  Percent_of_Ottawa_residents_80_1_dose
#>  Percent_of_Ottawa_residents_5_2_dose
#>  Percent_of_Ottawa_residents_12_2_dose
#>  Percent_of_Ottawa_residents_18_2_dose
#>  Percent_of_Ottawa_residents_30_2_dose
#>  Percent_of_Ottawa_residents_40_2_dose
#>  Percent_of_Ottawa_residents_50_2_dose
#>  Percent_of_Ottawa_residents_60_2_dose
#>  Percent_of_Ottawa_residents_70_2_dose
#>  Percent_of_Ottawa_residents_80_2_dose observed_census_ICU_p_acute_care
#> 
#> data:  VAR object var.a
#> Chi-squared = 68.027, df = 19, p-value = 1.951e-07